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Abstract 
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with Glauber dynamics. The scaling behavior of these quantities at low tempera- 
ture and large times is studied in detail. This analysis encompasses the self-similar 
domain-growth (aging) regime, the spatial and temporal Porod regimes, and the 
convergence toward equilibrium. The fluctuation-dissipation ratio admits a non- 
trivial limit value X^q = 1/2 at zero temperature, and more generally in the aging 
regime. 
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1 Introduction 



The ferromagnetic Ising chain with Glauber dynamics is one of the simplest non- 
equilibrium systems. Consider a finite system consisting of N Ising spins cr„ = ±1, with 
Hamiltonian 

= -J^cr„(j„+i. (1.1) 

n 

In its heat-bath formulation Glauber dynamics consists in picking, at every time step 

6t = 1/N, a site n = 1, . . . , at random, and updating its spin an{t) according to the 

stochastic rule i , j. u^oi, u\\ 

[+1 withprob. ^ + ^""^^^^^-^^^\ 

' ' l-tanh(/j/i„(t)) 
— 1 with prob. , 

where /3 = 1/T is the inverse temperature, and the local field hn(t) acting on (Tn(t) reads 

/in(t) = J(a„_l(t)+(T„+i(t)). (1.3) 

In the thermodynamic limit (A^ oo), each site is thus updated according to a Poisson 
process. 

At positive temperature, starting from a random initial condition, obtained for in- 
stance by a quench from infinitely high temperature, the system relaxes to its param- 
agnetic equilibrium state. The situation at zero temperature is different. Indeed, still 
starting from a random initial condition, the system is unable to relax to any of its two 
ferromagnetically ordered, symmetry-related, equilibrium states. On the contrary, do- 
mains of positive and negative magnetization grow forever, and, in the scaling regime, 
the system becomes statistically self-similar with only one characteristic length scale, 
the mean size of domains. This coarsening process is a consequence of the existence of 
spontaneous symmetry breaking 0, ||. 

The Ising chain is special because its critical temperature Tc is equal to zero. However 
its dynamical behavior illustrates that of generic statistical-mechanical models in the 
absence of quenched disorder. For the latter, starting from a random initial condition, the 
system relaxes exponentially to equilibrium in the high-temperature phase (T > Tc). At 
equilibrium, two-time quantities such as the correlation function C{t, s) or the response 
function R{t, s), where s (waiting time) is smaller that t (observation time), only depend 
on the time difference 

r = t-s>0, (1.4) 
and they are simply related to each other by the fluctuation-dissipation theorem: 

Re,iT) = -P^^. (1.5) 
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In the low-temperature phase (T < T^) the system undergoes phase ordering. In this non- 
equilibrium situation, C{t, s) and R{t, s) are non-trivial functions of both time variables, 
which only depend on their ratio at late times, i.e., in the self-similar scaling regime. 
This behavior is usually referred to as aging 0. Moreover, no such simple relation as 
eq. ( |1.5|) holds between correlation and response, i.e., R(t,s) and dC{t, s)/ds are no 
longer proportional. It is then natural to characterize the distance to equilibrium of an 
aging system by the so-called fluctuation-dissipation ratio ||, |[ 

ds 

In recent years, several works |Q, |^, |, 0, ||, ^ have been devoted to the study of 
the fluctuation-dissipation ratio. In the low-temperature phase of aging systems, such as 
glasses and spin glasses, or of systems exhibiting domain growth, X(t, s) turns out to be 
a non-trivial function of its two arguments. In the case of coarsening systems, analytical 
and numerical studies indicate that the limit fluctuation-dissipation ratio, 

= lim lim X(t,s), (1.7) 

vanishes throughout the low-temperature phase [0, |[. 

However, to date, only very little attention has been devoted to the fluctuation- dis- 
sipation ratio X{t,s) for non-equilibrium systems at criticality. From now on, we will 
only have in mind non-disordered systems. For instance one may wonder whether there 
exists, for a given model, a well-defined limit X^o at T = Tc, different from its value 
in the low-temperature phase. Indeed, a priori, for a system without disorder, such as 
a ferromagnet, quenched from infinitely high temperature to its critical point, the limit 
fluctuation-dissipation ratio X^o at T = Tc (if it exists) may take any value between 
Xoo = 1 (T > Tc. equilibrium) and Xoo = (T < T^. domain growth). 

The only cases of critical systems for which the fluctuation-dissipation ratio has been 
considered are, to our knowledge, the models of ref. |^ (random walk, free Gaussian field, 
and two-dimensional X-Y model at zero temperature) which share the limit fluctuation- 
dissipation ratio = 1/2, and the backgammon model, a mean-field model for which 
Tc = 0, where it has been shown that X^ = 1, up to a large logarithmic correction, for 
both energy fluctuations and density fluctuations [jlO|, |1T|. 

In the present work we investigate the non-equilibrium response function and fluc- 
tuation-dissipation ratio for the Glauber-Ising chain, both at zero temperature and in 
the scaling regime of low temperatures. Exact results for these two-time quantities are 
derived, exploiting the solvability of the model. We then perform a detailed analysis 
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of their scaling behavior at low temperature and large times. Since the computation 
of R{t,s) and of X{t,s) requires the knowledge of the equal-time correlation function, 
of the two-time correlation C{t,s), and of its derivative with respect to one of the two 



times [see eq. (|1.6| )], we will include the computation and the scaling analysis of the latter 
quantities, though a number of papers have already dealt with the study of correlations 
in the Glauber-Ising chain |T^, 0, |T3|, |TB|. Our presentation will thus be complete 



and self-contained, and the computation of all quantities will follow the same integral- 
transform methods. 

In a forthcoming publication |T7|, we will present a study of the critical response 
and fluctuation-dissipation ratio for the ferromagnetic spherical model in any dimension 
d > 2, and for the two-dimensional Ising model. Both models possess a whole low- 
temperature phase for T < Tc, hence they are faithful representatives of generic domain- 
growth systems. 

One salient outcome of these joint works is the realization that the limit fluctua- 
tion-dissipation ratio is a novel universal characteristic of critical dynamics (see the 
discussion at the end of the present paper). For the Ising chain at its critical temperature 
Tc = 0, we obtain the value Xoo = 1/2. The occurrence of a common value for Xoo 
between the Ising chain and the models considered in ref. seems rather coincidental, 
for the time being. 

In the following, we shall make use of the parameters fi and 7, defined by 

u = -lntanh(/3J), 7 = tanh(2/3J) = — (1.8) 

cosh/x 



These parameters are related to the equilibrium correlation length ^eq |@| and to the 
equilibrium relaxation time Tgq [see eq. ( f^.l6D ] by 

1 1 

Ceq = ; Teq = :; • (1-9) 

/i 1-7 

Both diverge exponentially fast as the critical temperature Tc = is approached, accord- 
ing to 

U - ^eq - ^ - (1.10) 

hence we have the scaling law 

req ^ (1.11) 

corresponding to a dynamical critical exponent 2; = 2, and reflecting the diffusive nature 
of domain growth at zero temperature. 

We close up this introduction with a discussion about relevant time scales. For two- 
time quantities such as C(t, s), -R(t, s), X{t^ s), three time scales are to be compared in 
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the scaling regime of long times and low temperature, namely s, r = t — s, and Teq. 
Considering any two of these time scales as comparable, with an arbitrary ratio between 
them, and small (or large) compared to the third one, six different regimes can be defined 
a priori. Three of them will be of interest in this work, which we summarize here for 
convenience: 

T ~ Tcq -C s : Equilibrium (1.12) 

s ~ r ^ Tcq : Self-similar domain growth (aging) (1-13) 

T <^ s ~ Teq : Early-time or temporal Porod regime (1-14) 

Aging persists forever at zero temperature, just as the coarsening process itself, while it 
is interrupted at the time scale Teq at any low but finite temperature. 



2 Magnetization profile 

In this section we present the methods used throughout the paper on the example of the 
relaxation of the magnetization profile 

M„(t) = (cr„(t)) , (2.1) 

where the brackets denote an average over the thermal history of the system and over 
the initial conditions, which are unspecified for the time being. 

The time evolution of the magnetization M„(t) is readily deduced from the stochastic 
rule ( pr2|) , and is given by 

^Mn{t) _ ^ {(3K{t))) . (2.2) 



dt \ 

Since the local field hn only assumes three symmetric values, and ±2J, we have 

tanh(/?/i„) = — hn= 2 (^"-i + ^"+1)' (2-3) 
leading to the linear evolution equations 

'^MM = -Mr,{t) + |(M„_i(t) + M„+i(t)). (2.4) 

The existence of closed, linear equations for the magnetization, and more generally for 
higher-order correlation functions, ensures the solvability of the model |jl[. 

In order to solve coupled equations of the form ( p.4|) , we shall make an extensive use 
of Laplace and Fourier transforms. For any quantity fn{t), depending both on continuous 
time t and on the discrete site label n, we introduce 
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the temporal Laplace transform 

fnip) = /f m e-^* dt, m = / ^ e^ 
the spatial Fourier transform 



(2.5) 



-mq 



fn{t)=l ^f{q,t)e'^^, 



27r 



(2.6) 



the Fourier-Laplace transform 
^ Jo 



Using the above integral transforms, eq. (|2.4|) can be recast as 

dMP(g,t) 



dt 



(7cosg-l)M^(g,t), 



or else as 



yielding 



pM^^{q,p) = {-f cos q- I) M^^{q,p) + M^{q,t = 0), 

M^{q,t = 0) 



M^L(g,p) 



(2.8) 
(2.9) 
(2.10) 



p + 1 — 7 cos q 

Consider first the locally magnetized initial condition where the spin at the origin is 
pointing upwards, i.e., ao(t = 0) = +1, but the configuration is otherwise totally random. 
The corresponding solution Gn{t) to eq. ( |2.4] ) is the Green's function of the problem. We 
have Gn{t = 0) = 5„,o, so that G^{q,t = 0) = 1. Eq. {^J^ then reads 



G'^iq,p) 



1 



hence 



and 



p + 1 — 7 cos q 



Gn{t) = e- 



dq 
2^ 



^'yt cos q+inq 



(2.11) 

(2.12) 
(2.13) 



where is the modified Bessel function. Eq. ( p. Ill) also yields 



Gl:ip) 



(p+l)2-72 



p + l-^(p+l)2-72 

7 



\n\ 



(2.14) 
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Coming back to eq. (p^), its general solution in direct space is obtained by inverting 
eq. (|2.10|) , yielding the following spatial convolution 



M„(t) = M„(t = 0) * Gn{t) = J2 Mm{t = 0)Gn~m{t)- (2.15) 

m 

In particular, for a uniformly magnetized system, i.e., M„(t = 0) = M for all n, we have 
an exact exponential relaxation 

^ = E G„(t) = G^iq = 0,t)= e-/-^ (2.16) 

n 

at any finite temperature, where the relaxation time Teq is given by eq. ( |1.9| ). 

Throughout the following, we shall be mostly interested in the non-equilibrium scaling 
regime of long times, large distances, and low temperatures, such that t and Teq are 
simultaneously large but comparable, i.e., their ratio t/xeq is arbitrary, and the same 
holds for n and ^eq: 

t ~ Teq > 1, n ~ ^eq > 1- (2.17) 

The scaling law ( |1 . 1 1| ) then implies 

t ~ Teq ~ ^2 ~ ^eq- (2-18) 

As a consequence, p ~ ~ /x^ are simultaneously small, and eq. (p.ll| ) simplifies to 



G^'^iq^p)^-—^^—^. (2.19) 
2p + + fi^ 



Performing the inverse Laplace transform first, we obtain 



t n^^ 



G,.M«-^exp^----j. (2.20) 



which is the scaling form of the Green's function (|2.13|) , involving the variables t/req and 
n/ \/t, in agreement with eq. ( |2.18| ). At zero temperature the exponential damping factor 
exp(— t/req) is absent, so that Gn{t) is a function of n/\^ only, reflecting the statistical 
self-similarity of the coarsening process. In this case the Gaussian profile is normalized, 
emphasizing again the underlying diffusive mechanism. 



3 Equal-time correlation function 

The equal-time correlation function between any two spins is defined as 

C(m,n,t) = (a„(t)(T„(t)). (3.1) 
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In the following, unless otherwise specified, we shall consider a random initial condi- 
tion, defined by averaging over all possible initial configurations with equal weights. This 
procedure corresponds to quenching the system at time t = from its equilibrium state 
at infinitely high temperature. Because the invariance under spatial translations along 
the chain is preserved by the dynamics, the correlation function C{m,n,t) only depends 
on the distance |n — m| between both spins. We denote it as 

C(m,n,t) =C„_„(t). (3.2) 

We have in particular 

Co{t) = 1. (3.3) 



The equal-time correlation function can be shown, in analogy with eq. (|2.4| ), to obey 
coupled linear differential equations of the form [|T| 

^Ml = _2C7„(t) + 7(a-i(t) + a+i(t)) (n ^ 0), (3.4) 
at 

with the condition ( p.3|) , and the initial value C„(t = 0) = 5„_o- 

In order to solve eq. ( p.4| ), we complete them by the corresponding equation for n = 0, 
with a time-dependent source v{t) in the right-hand side, to be determined in such a way 
that the condition (|3.3| ) be fulfilled. In other words we consider the equations 

= -2Cr.it) + 7(C„_i(t) + C„+i(t)) + v{t)6,,,o, (3.5) 

at 

which read, in Fourier-Laplace space, 

p(:7^L(g,p) = 2(7cosg-l)C^L(g,j9)+t;L(p) + l, (3.6) 



yielding 

p + 2 — 27 cos q 



c-fe,p) = :y:' . (3.7) 



The condition (|3.3|) then reads 

r^C--iq,p) = ^£M±^ = \ (3.8) 
Jo 271 ^(p + 2)2_4y p 



hence 



We thus obtain the result 



'(p + 2)2 -472 

v^{p) = ^ - 1. (3.9) 

p 



.PL, . ,/(p + 2)2-472 



p{p + 2 — 27 cos q) 
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or else 

C'^iP) = -\^- ^ (3-11) 

which we now discuss. 
• Equilibrium 

At any finite temperature, and for long enough times {t ^ Teq), the correlation func- 
tion Cnit) converges to its equilibrium value, C„,eq = limp_>o(pC^(p)). From eq. 
we recover the well-known expression 



• N on- equilibrium 

Let us analyze the non-equilibrium properties in the regime ( |2.18|) of low temperatures 
and long times. In this regime the expressions ( ^.10 ) and ( 3.1 1| ) scale as 



C'\^^P)^^^^W^, (3.13) 
and 

^ • (3.14) 

P 

Eq. (A.l) leads to the following explicit scaling form 

Cnit) ^ \ je^l^l erfc(^^ + ^v^^ + e"'^!"! erfc(^^ - ^iVt^ | (3.15) 

for the equal-time correlation function, involving two scaling variables, /x^t ~ 2t/req and 
n/y/t. The error function erf z and the complementary error function erfc^; are defined 
as 

erf z = —;= e ^ dx, eifcz = 1 — erf 2; = —= / e dx. (3.16) 

Jo y/TT Jz 



An alternative expression for Cnit) is 

which can be obtained either from eq. ( |3.13| ) or from eq. ( |3.15| ). The following limiting 
situations are of interest. 

(i) At zero temperature, and more generally in the domain-growth regime (1 <^ t <^ Teq), 
we can set /x = in eq. (|3.15| ), which thus simplifies to 



Cnit)^ erfc^f^ . (3.18) 



9 



This expression |T5| involves only one scaling variable, n/\/t, reflecting the fact that, 
in this regime, the pattern formed by the domains is self-similar. 

(ii) In the opposite regime {t ^ rgq), the expression p.l5| ) converges exponentially fast 
to its equihbrium value ( p.l2| ), as 

I I / , \ -3/2 

C.(t)-C.,eq^--^ — e-2*/-''. (3.19) 
4v27r 

(iii) The short-distance or Porod regime 

K |n| < ^eq ~ Vt (3.20) 
is yet another situation of interest 0. In this regime, eq. ( p.l5|) becomes 



Cn{t)^l-A(t)\nl (3.21) 

with 

A{t) = ^= + /i erf(/iv^). (3.22) 



This result may also be found from eq. ( |3.14|) , by noticing that |n|\/p+~/? ^ 1, and 
using eq. (A. 2). The expression ( |3.22|) for the amplitude A{t) interpolates between the 



power-law behavior 

A{t) ^ ^7 (3-23) 



in the domain-growth regime, in agreement with eq. (|3.18| ), and the equilibrium value 



Aq ~ /i, (3.24) 

in agreement with eq. ( ^.12 ). 



The result ( |3.21|) also holds for finite values \n\ = 0,1,2,... of the distance between 
spins, in the scaling regime (,^eq ^ 1, t ^ 1). In particular the density of defects (domain 
walls) in the system reads 

. « m. (3.25) 



The result (|3.21| ) is thus in quantitative agreement with the general predictions given in 



ref. on the so-called Porod singularities of correlation functions, whose form depends 
on the nature of topological defects, and whose amplitude is known in terms of the density 
Pdef of these defects. 
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The spatial range over which ferromagnetic order has propagated at time t can be 
measured by the dimensionless susceptibihty 

xW = ECn(t) = C7^(g = o,t), (3.26) 

n 

for which eq. ( p.lO| ) yields 

. '-J^I±S, (3.27) 
Throughout the non-equilibrium scaling regime ( |2.18|) , we have 

x'-ip) ^ r^^ , (3.28) 

i.e., using eq. (A. 3), 

X(t) ^ - erf LVt) . (3.29) 
This result interpolates between the square-root behavior 

X{t) ^ 4^ (3.30) 
in the domain-growth regime {t <^ Teq), in agreement with eq. (|3.18|) , and the limit 

Xeq ~ - (3.31) 

as equilibrium is approached (t S> Xcq). This last expression is the scaling behavior of 
the exact equilibrium value 



= cotanh t = " + = e^^\ (3.32) 

2 7 - 1 + Vl - 7^ 

corresponding to eq. ( p.l2| ). 

Finally, the dimensionless product 

g(t) = A{t)x{t) = 2 erf(^Vt) + erf (/iVt) (3.33) 

is a form factor characteristic of the correlation profile. It increases from the value 



Q = A/tc in the domain-growth (aging) regime [see eqs. ( p.23| ), (|3.3CI|) ] to the equilibrium 
value Qcq = 2 [see eqs. j ^M ), i^M)]- 
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4 Two-time correlation function 



We now consider the two-time correlation function 

C(m, n, t, s) = (a„(t)(Tm(s)) , (4.1) 

where s is the waiting time and t = s + r is the observation time. For a random initial 
condition, invariance under spatial translations yields 

C(m,n,t,s) = C„_„(t,s). (4.2) 



This two-time correlation function can be shown, again in analogy with eq. ( |2.4|) , to 
obey the coupled linear partial differential equations 

^^^I^ = -C4t, s) + |(a_i(t, s) + C„,+i(t, s)) (4.3) 
for t > s, with the initial value 

C„(s,s) = a(s), (4.4) 

at time t = s, i.e., r = 0, where the right-hand side is the equal-time correlation function, 
given in Fourier-Laplace space by eqs. ( |3.10| ) and ( ^.11| ). 



The second argument s plays the role of a parameter in eq. ( |4.3| ), hence this equation 
is formally identical to eq. (|2.4|) , with the initial condition ( [4.4|) playing the role of 



Mn(t = 0). The solution of eq. (|4.3|), seen as an evolution equation in the r variable, is 
therefore formally identical to that of eq. (p.4|) , and reads 



Cn{s + T,s) = Cn{s) * G„(r) = C„ (s)G„_^(r) . (4.5) 

m 

In order to write down this solution explicitly, it is convenient to introduce the double 
Laplace transform of the function C„(s + r, s), where ps is conjugate to s, the transform 
being denoted by L^, and p is conjugate to r, the transform being denoted by L: 

roc roo 

C]i'^'{p,Ps)= / C„(s + r,s)e-(^'=^+P^)dsdr. (4.6) 

JQ Jo 

With this definition, the solution of eq. (^^ ) reads 



g"""-(g,P,P.)= f""^^'^-^ , (4.7) 

p + 1 — 7 cos q 



e., using eq. ( ^.lOj ) 



C'^^^{q,p,Ps) = . % ^ r. (4.8) 

Ps{Ps + 2 — 27 cos g) (p + 1 — 7 cos q) 



12 



In the following, we shall mostly be interested in the diagonal component of the cor- 
relation function, or autocorrelation function, Co(t, s). Eq. (|4.8| ) yields, upon integration 
over q, 

1 



PsiPs-'^P)\\ (p+l)2-72 



(P' + ^r-^^'-'.2\. (4.9) 



The discussion given at the end of section 3 on the behavior of C„(t) in various regimes 
is readily extended to the two-time correlation function C„(t, s) considered in the present 
section, using its integral-transform expressions (^^ ), (|4.9| ). 

• Equilibrium 

The equilibrium correlation function is obtained by letting s — oo (in practice, s 3> 
Teq), while keeping r = t — s fixed. Eq. (|4.8| ) simplifies to 



Ceq = 7^ -y (4.10) 

^ (1 — 7 cosg)(p + 1 — 7 cosgj 

which yields, performing the inverse Laplace transform first, 

C„,eq(r) = e-^l"l - ^^-1^ [ Gn{u) du = v/l-7' Gn{u) du. (4.11) 

In particular the equilibrium autocorrelation function reads 

Co,eq(r) = ^1-7' r G'o(m) du. (4.12) 



In the scaling regime, where r and Tgq are large and comparable, using eq. (|2.2CI| ), the 
above formula simplifies to 

Co,eq(r)^ erfc(^^^. (4.13) 

• N on- equilibrium 

The generalization of the non-equilibrium scaling regime (|2.18|) of long times, large 
distances, and low temperatures, to the present case of two temporal variables consists in 
taking s, r, and Teq simultaneously large but comparable, with arbitrary ratios between 
them, and the same for n and i^eq, i-e., using again eq. (|1.11|) , 



s ~ r ~ Teq ~ n ~ Q > 1. (4.14) 

As a consequence, ps ~ p ~ ~ /i^ are simultaneously small. In this regime, eq. ( |4.9| ) 
reads 

C^^'{p,Ps) ~ „ , , - — , (4.15) 

PsVWTJ? [VpTTJ^ + V^pTJ?) 
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Performing first the inverse Laplace transform of this expression with respect to p, using 
eq. (A. 4), we get 

C^^iT,ps) ^ erfc (^ips + i^^)T/2) . (4.16) 

Performing then the inverse Laplace transform with respect to ps, using eq. (A. 5), we 
obtain 



\/2t 9 , /■« e~^ " 
Co{s + T,s)^^e-'^'^/' du. (4.17) 

TC Jo [2u + T)y/u 



The change of variable u = (t/2) tan^ 6 yields 

2 /•© 



Co{s + T,s)^- f e-'^'"/(2co.-e) ^^^^^ 

TT JO 



with 

e = arctanW— . (4.19) 
Eq. ( 4.18| ) gives the general scaling form of the two-time autocorrelation function. It can 



be alternatively expressed as a double Gaussian integral 



4 r'^ /•5-v/2s/t „ „ 

Co(s + r,s)^-/ e-^^+^Mry. (4.20) 



TT Jny/T/2 JO 

Only in the symmetric situation O = 7r/4, i.e., r = 2s, or t = 3s, this integral simplifies 



to 



Co(3s,s)^^(l-erf2(/iv/j)). (4.21) 



Yet another reformulation of the result ( [4.18|) reads 



where $i is the confluent hypergeometric series in two variables [jl9[, namely the scaling 
variables 2s/r and /i^s ~ 2s/req. 



The fluctuation-dissipation ratio ( |1.6| ) involves the derivative of Co(t, s) with respect 
to s at fixed observation time t, for which eq. (|4.18 ) yields the simpler expression 



9Co(t = s + r,s) /2(s + r)e-^'^ ^ \ e-'^Va 

^ ' ^ / V ; +^erf(/iv^) ^^=. (4.23) 



ds y 2s + r a/tts y v27rr 

Let us discuss some limiting cases of eqs. ( |4.18| ) and ( [4.23| ). 



(i) At zero temperature, and more generally in the domain-growth regime ( |1.13|) , the 
autocorrelation function exhibits aging, i.e., it only depends on the ratio t/s. Eq. ( [4.18| ) 
indeed simplifies for = to Co(t, s) ~ 26/7r, i.e., 

2 [2s 

Co(s + r, s) ^ -arctanW — . (4.24) 

TT V "7" 
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This expression O, |T^, |T5|, |T6i assumes the scahng form 



with 



and 



Co{s + T,s) ^ F{x) 



X=- = 1 + ->1, 

s s 



^/ N 2 2 

F[x] = — arctan W . 

vr V a; — 1 



We have F{3) = 1/2, in agreement with eq. (|4.21| ). Similarly, eq. (|4.23| ) becomes 

dCo{t,s) 



ds 



(4.25) 
(4.26) 
(4.27) 

(4.28) 



with Fi{x) = —xdF{x)/dx, i.e.. 



Fi(x) 



X 



tt{x + 1) V a; — 1 



(4.29) 



(ii) For s oo with r fixed, i.e., Ps — *• 0, eq. ( [4.18|) converges to the equilibrium 
form ( |4.13| ); this is obvious using eq. ( ]4.16| ). 

(iii) The early-time regime, or temporal Porod regime (|1.14|) , is the counterpart of the 
spatial Porod regime i WM - From eq. (14231 ), we find 

dCo{s + T,s) A{s) 



ds VSttt' 

where A{s) has been defined in eq. ( |3.22|) . We thus obtain, to leading order, 

[2r 

Co{s + T,s) ^ l-A{s) 



(4.30) 



(4.31) 



The replacement of \n\ in eq. ( p.21| ) by (a constant times) ^/T in eq. ( |4.31| ) refiects once 
more the underlying diffusive mechanism. Finally, the behavior (|3.23| ) of A{s) in the 
domain-growth regime is in agreement with eq. ( [4.24| ), while eq. ( p. 24 ) matches eq. ( |4.13| ). 
The generalization to the spatio-temporal Porod regime 



1 n ~r^s~r. 



eq 



(4.32) 



is straightforward. One has, using eqs. ( [4.5|) , (t3.21|) , (|2.20|) 



Cnis + T,s)^il-A{s)\n\)* 



-n2/(2r) 



V2 



TCT 



(4.33) 
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where the convolution involves the spatial coordinate n. Replacing the discrete sum by 
an integral, we obtain 



Cnis + T,S)^1- A{S) { e-"'/(2r) + ^ g^f 



(4.34) 



This result interpolates between the spatial behavior (|3.21|) of the equal-time correlation 
function, for r = 0, and the temporal behavior ( [4.31| ) of the autocorrelation function, for 
n = 0. 



5 Two-time response function 

Suppose now that the system is subjected to a small magnetic field if„(t), depending on 
the site label n and on time t > in an arbitrary fashion. This amounts to adding to 
the ferromagnetic Hamiltonian ( |1.1|) a time-dependent perturbation of the form 

Sn{t) = -J2Hn{t)an{t). (5.1) 

n 

The dynamics of the model is still given by the stochastic rule ( |1.2|) , where the local field 
hn{t) acting on the spin cr„(t) now reads 

Kit) = J{an-i{t) + an+i{t)) + Hn{t). (5.2) 



We again consider a random initial condition. Causality and invariance under spatial 
translations imply that the magnetization M„(t) at time t reads 

Mn{t) = ((J„(t)) =(3 f duY, Rn-m{t, u)Hm{u) + ■ ■ • , (5.3) 

to first order in the magnetic fields Hn{t). This formula defines the two-time dimension- 
less response function Rn-m{t,s) of the model. A more formal definition reads 



SMn{t) 



6Hrr,(s) 



(5.4) 

{H=0} 



The evolution equation ( |2.2| ) still holds in the presence of arbitrary external magnetic 
fields. Furthermore we have, to first order in the magnetic field 

tanh(^/i„) = tanh(/3J((J„_i + o-„+i)) + (3Hn (l - tanh^(/3J(cr„_i + cr„+i))) H , (5.5) 
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together with the identities 



tanh(/3J((j„_i + an+i)) = l^i^^n-i + c^n.+i), 

tanh^(/5J(o-n-i + o-„+i)) = y (1 + a„_ia-„+i). (5.6) 



Inserting these expressions into eq. (|2.2|) , we readily obtain that, again to first order 
in the magnetic field Hn{t), the magnetizations Mn{t) obey inhomogeneous differential 
equations of the form 

= -M„(t) + |(M„_i(t) + M„+i(t)) + (^1 - ^(1 + + • • • 

(5.7) 

As a consequence, the two-time response function Rn(t,s) itself obeys coupled linear 
differential equations, of the form 

dR it s) ^ 

^ = -i?n(t, s) + ^(i?„-i(t, s) + i?„+i(t, s)), (5.8) 

for t > s, with the initial value 

Rn{s, s) = w{s)6n,o, (5.9) 

and with 

w{s) = l-^{l + C2is)), (5.10) 

where 6*2(5) is the equal-time correlation two sites apart, which is given in Laplace space 
by eq. ( p.ll| ). Eq. (|578|), with its initial condition ( |5.9| ), is formally identical to eq. ( |4.3| ) 
with initial condition ( |4.4| ). Hence its solution reads 

Rn{s + T,s) =w{s)Gn{T), (5.11) 

or, in Laplace space with respect to s, 

i?^^(r,p,) =w;L^(p,)G„(r). (5.12) 



Using eq. (|3.11|) , we obtain 



w'^'iPs) = ^^(^. + 2)2-47^ - 1 - f • (5.13) 

Let us now discuss the general expression of the response function, given in Laplace 
space by eqs. ( |5.12|) and ( ^.13|) , along the same lines as for the two-time correlation 
function. 
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• Equilibrium 

At equilibrium, i.e., for s ^ oo with r fixed, one gets, either using eq. ( p.l2| ) or taking 
the — > hmit of eq. ( |5.1ci| ), 



«;eq=Vl-7^ (5.14) 

hence 

i?n,eq(r) = ^/l^Gnir). (5.15) 

Comparing this equation to eq. (|4.11| ) leads to the identity 

Rn,e^{T) = (5.16) 



which is the fiuctuation-dissipation theorem ( |1.5| ), in dimensionless form. 
• Non- equilibrium 

For simphcity we again focus our attention on the diagonal component Ro(t,s). In 
the scaling regime ( [4.14| ), we have 

w '{Ps) = • (5.17) 

Ps 

Eq. (A. 2) implies that the function w{s) coincides in this regime with the amplitude 
A{s), defined in eq. (p.22| ). Using eq. ( |5.12|) , and the scaling form ( |2.20| ) of Go{t), we 
obtain the remarkably simple result 

Rois + T,s)^ ' . (5.18) 

This general scaling form of the two-time response function further simplifies in the 
following limiting cases. 

(i) At zero temperature, and more generally in the aging regime ( |1.13| ), the result ( p. 18] ) 
becomes 

R,{s + T,s)^^=. (5.19) 
7rv2sr 

We thus obtain a scaling form similar to eq. ([4.28|) , namely 

R,(s + T,s)^^, (5.20) 

with 

F2{x) = . ^ (5.21) 



7r^2(a;-l) 

(ii) In the early-time regime ( |1.14|) , we obtain 



V27rr 

where the right-hand side is identical to eq. (|4.3CI| ). 



i?o(s + r,s)^^, (5.22) 
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6 Fluctuation-dissipation ratio 



As mentioned in the introduction, a way of characterizing the violation of the fluctua- 
tion-dissipation theorem (|1.5|) in a non-equilibrium situation consists in introducing the 
fluctuation-dissipation ratio X(t, s) of eq. (pTBI). In the present situation, we set 

X(t,.) = -?#^, (6.1) 
^' ^ dCo{t,s) ^ ^ 

ds 

because the response function Ro(t, s) is dimensionless. 

In the non-equilibrium scaling regime (CT) , eqs. (14231) and (|l8|) yield 

X{s + T,s)^ , e-^-- + /.Aierf(;zVi) ^g^) 

2(s + r) 2, , ^ ci r\ 
■e~/^ * -|- ii^T^s ert(/XY's) 



2s + r 

which is the scaling form of the fluctuation-dissipation ratio in the variables s/rgq and 
r/Tcq. We again discuss some limiting cases. 

(i) For s Teq, with r and Teq fixed, the fluctuation-dissipation ratio converges exponen- 
tially fast to its equilibrium value X^^ = 1, according to 



(ii) The result (|6.2| ) also shows that we have X{t, s) ~ 1 in the early-time regime ( |1.14| ), 
consistently with the identity between expressions ( [4.30|) and (|5.22| ). Notice, however, 
that the initial value X{s,s) of the fluctuation-dissipation ratio is not identically equal 
to 1. For instance, at zero temperature, eqs. ([4.9|) and ( |5.13|) yield after some algebra 

2s(/o(2s) +/i(2s)) 

This expression behaves as X{s,s) ~ 1/2 in the regime s <^ 1 of no physical interest, 
while it converges to unity after a microscopic transient regime, as X{s, s) ~ 1 — l/(4s) 
for s ^ 1. 

(iii) At zero temperature, and more generally in the aging regime (|1.13|) , eq. ( |6.2|) sim- 
phfies to 

X{s + T,s) ^ A:{x), (6.5) 
where the scaling function X{x) reads 

Mix) = (6.6) 
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This result is consistent with the scahng laws (^4.28| ) and ( p. 20 ), as we have 



X x) 



FAx) 



(6.7) 



The fluctuation-dissipation ratio decreases from the initial value ^(0) = 1, corresponding 
to equilibrium, to the non-trivial asymptotic value 



= X{oo) = 1/2. 



(6.8) 



Recent developments on aging systems 0] suggest the following alternative presenta- 
tions of the above results concerning the aging regime. 

Firstly, there is a functional relationship X{C) between the fluctuation-dissipation 
ratio X = X{t, s) and the two-time correlation function C = Co{t, s) in the aging regime, 
as a consequence of eqs. ( [4.25|) and (|6.5|) . Indeed, eliminating the time ratio x between 
eqs. ( [4. 271 ) and (|6.6|), we obtain 

1 



XiC) 



TX 



COS^ 



(6.9) 




Figure 1: Plot of the relationship X[C) of eq. (6^) between fluctuation-dissipation ratio 
X(t, s) and correlation function Co(t, s) in the aging regime. 



Secondly, the dimensionless integrated response function 



Po(i,s) 
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Roit, u) dw. 



(6.10) 



proportional to the thermoremanent magnetization M, reads 



Po{t,s) = / X{t,u)dCo{t,u) 



(6.11) 



using the definition ( |6.1| ) of X{t,s). In the aging regime, the existence of the functional 
relationship (|6.9| ) implies 

Po{t,s)^piCoit,s)), (6.12) 



with 



i.e., explicitly. 



p{c)= r x{c')dc\ 

Jo 



p[Cj = arctan I ^an ] ] ■ 



(6.13) 
(6.14) 



1.0 



0.8 - 



0.6 




0.4 



0.2 - 



Figure 2: Plot of the relationship p(C) of eq. (|6.14|) between integrated response po{t, s) 
and correlation function Co(t, s) in the aging regime. 



The functions X{C) and p{C) are respectively shown in Figures 1 and 2. For C ^ 1, 
i.e., X = t/s 1, we have 



TT 



XiC) = l--{l-C)' + 



p{C) = ^-{l-C) + ^{l-C)' + ---, (6.15) 



while for C ^ 0, i.e., x — > oo, we have 
^ ^ 2 16 



(6.16) 
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Let us point out that the slope of the p{C) curve near the origin is given by the asymptotic 
value X{C = = = 1/2. 

It is clear from eq. ( |6.2| ) that the fluctuation-dissipation theorem is maximally violated 
for r ^ s, at least in the scaling regime (|2.18| ). In this situation, eq. ( |6.2|) yields the 
prediction 



for the asymptotic fluctuation-dissipation ratio 



Xjs] 



lim X(s + T, s) 

T^OO 



(6.17) 



(6.18) 




0.2 



0.5 



S/T 



1.0 



eq 



1.5 



Figure 3: Plot of the asymptotic fluctuation-dissipation ratio Xas(s) in the non-equili- 
brium scaling regime, as given by eq. ( |6.17| ), against s/r^q. 



Figure 3 shows a plot of the prediction ( |6.17| ) for Xas(s), against the ratio s/req. In 
the aging regime (1 ^ s ^ Tcq), the asymptotic fluctuation-dissipation ratio smoothly 
departs from its limit value X^o = 1/2, as 



XJs) ^ - + — + 



(6.19) 



eq 



while it converges exponentially fast to its equilibrium value Xeq = 1 for s ^ Teq. 

X,,(s) ^1 ,, 

2ns 



(6.20) 
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7 Discussion 



This work is devoted to the non-equihbrium dynamics of the ferromagnetic Ising chain, 
quenched from infinite temperature to finite temperature, and evolving under Glauber 
dynamics. We have exploited the solvability of this model in order to derive exact ex- 
pressions for the spin autocorrelation function Co(t, s) = {(To(t)(To{s)), and the associated 
response function RQ{t,s) and fluctuation-dissipation ratio X{t,s), with t (observation 
time) > s (preparation time) and t = t — s. While our study of correlations complements 
a well-investigated field, the results concerning the response and the fluctuation-dissipa- 
tion ratio are entirely novel. 

We have given a detailed analysis of the scaling regime of low temperatures, large 
distances (proportional to the equilibrium correlation length ^eq) and large times (pro- 
portional to the relaxation time Teq ~ ^gq). This scaling regime encompasses several 
limiting situations of interest: 

• Self-similar domain- growth or aging regime (1 ^ s ~ r ^ ^cq)- At zero temperature, 
because of the self-similarity of the coarsening phenomenon, the various two-time observ- 
ables obey simple scaling laws involving only the time ratio x = t/s. This situation is 
typical of an aging system [^. The corresponding scaling functions for the autocorrela- 
tion function Co(t, s), its derivative dCo(t, s)/ds, the response function -Ro(^; s), and the 
associated fluctuation-dissipation ratio X(t,s) are given explicitly in eqs. ( [4.27| ), ( |4.29| ), 



( ^.21| ), and ( |6.6|) . These functions are universal, implying in particular that they are un- 



changed if the initial condition contains short-range correlations. These results imply the 
existence of a non-trivial relationship X{C) [see eq. ( |6.9| )] throughout the aging regime. 
At any low but finite temperature, the self-similar domain growth and the associated 
aging phenomena are interrupted for times of order Xeq. 

• Spatial (1 ^ \n\ ^ ~ y/s) and temporal (r ^ s ~ Teq) Porod regimes. In these 
regimes, the autocorrelation function departs from its value of unity in a singular fashion, 
involving either an \n\ [§] or a ^/T dependence [see eqs. (|3.21|) , (|4.31| )]. We have derived 



a prediction p.22| ) for the common prefactor A{t) of both these laws, as well as the 
interpolation formula ( f4.34|) in the spatio-temporal regime. 

• Equilibrium regime (r ~ Tgq <^ s). At low but finite temperature, the system converges 
toward equilibrium for times larger than Teq. The exponential law of convergence has also 
been studied in detail [see eqs. ( p.l9|) , ( |6.3|) , (|6.20|) ]. 

As already underlined in the introduction, one of the most salient outcomes of this 
study is the non-trivial limit value X^o = 1/2 of the fluctuation-dissipation ratio for 
1 ^ s ^ t in the aging regime. The value of X^o for a system quenched to its critical 
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point could be a priori any number between = 1 (T > T^. equilibrium) and X^ = 
(T < Tc'. domain growth). It turns out that the answer is exactly half-way between these 
bounds for the Glauber-Ising chain. A forthcoming work [|T7| is devoted to the limit 
fluctuation-dissipation ratio in systems having a finite T^, and exhibiting domain growth 
in a whole low-temperature phase. The spherical model in any dimension d > 2 and the 
two-dimensional Ising model will serve as benchmarks and lead us to claim that X^o is a 
novel universal characteristic of critical dynamics, intrinsically related to non-equilibrium 
phenomena, which can take any value, at least in the range < X^o < 1/2. 

Let us anticipate that the limit fluctuation-dissipation ratio X^o appears as an am- 



plitude ratio in the sense used in critical phenomena. For a critical quench in such 



generic models, we have indeed 



Cit,s) ^ s-(''-'+'')/^F(x), 
T,R{t,s) ^ s-^-('^-2+W-^2(x). (7.1) 



These scaling laws generalize eqs. ( [4.25| ), ( [4.28| ), ( ^.20| ). The dynamical exponent z is 



accessible from the study of dynamical critical phenomena at equilibrium, while x = t/s 
is again the time ratio. The fluctuation-dissipation ratio therefore scales as 

Xit,s)^X{x) = ^^, (7.2) 

where the scaling function X{x) is again universal. For large values of this ratio {x ^ 1, 



i.e., 1 <^ s <^ t), the scaling functions entering eqs. ( [7. 1|) fall off as 



F 


[x) 


^ Bx- 


Xc/Z 


Fi 


[x) 


^ Bix 


-Xc/z 


F2 


[x) 


^ B2X 


~Xc/z 



(7.3) 

where Ac, related to the critical initial-slip exponent 6c |^ by Ac = d — z6c, is a dynamical 
critical exponent which only manifests itself in non-equilibrium phenomena The 
exact results derived in this paper for the Ising chain agree with the above scaling laws 
for = 1, z = 2, Ac = 1, and 9c = 0. The limit fluctuation-dissipation ratio 

X^ = X{oo) = ^ (7.4) 

thus appears as a dimensionless amplitude ratio associated with the behavior ( [7.31 ). It is 
therefore a novel universal quantity in non-equilibrium critical dynamics. 
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A Some useful inverse Laplace transforms 



In this Appendix we list a few inverse Laplace transforms, which can be found in ref. 



and have been used in the body of this article. Notations are as in eq. (|2.5| ). The symbols 
a, b denote complex parameters with a positive real part. 



/"(P) fit) 



P 
a 



p\/p + a? 



1 



e""^ erfcf + aM + e""'' erfcf - aM (A.l) 



P \2Vt J \2Vt 



T) I Qj^ pt ^ ^ 

—= + a erf (a^Ti) (A.2) 



erf (a^Tt) (A.3) 



erfc (yap) (^-5) 
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